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ABSTRACT 

High spatial resolution observation of an emerging flux region (EFR) was 
done using a vector magnetograph and a Ha Lyot filtergraph with Domeless 
Solar Telescope at Hida Observatory on October 22, 2006. In Ha wing images, 
we could see many Ellerman bombs (EBs) in the EFR. Two observation modes, 
slit scan and slit fixed, were performed with the vector magnetograph, along 
with Ha filtergraph. Using Ha wing images, we detected 12 EBs during slit 
scan observation period and 9 EBs during slit fixed observation period. With 
slit scan observation, we found that all the EBs were distributed in the area 
where the spatial gradient of vertical field intensity was large, which indicates 
the possibility of rapid topological change in magnetic field in the area of EBs. 
With slit fixed observation, we found that EBs were distributed in the areas of 
undulatory magnetic fields, both in vertical and horizontal components. This 
paper is the first to report the undulatory pattern in horizontal components 
of magnetic field, which is also evidence for emerging magnetic flux by Parker 
instability. These results allow us to confirm the association between EBs and 
emerging flux tubes. Three triggering mechanism of EBs is discussed with respect 
to emerging flux tubes: 9 out of 21 EBs occurred at the footpoints of emerging 
flux tubes, 8 out of 21 EBs occurred at the top of emerging flux tubes, and 4 out 
of 21 EBs occurred at unipolar region. Each case can be explained by magnetic 
reconnection in the low chromosphere. 
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INTRODUCTION 



Ellerman bombs (EBs) are short-lived, small-scale bright points observed best in the 
wings of chromosp heric lines, such as Ha and Ca II H line. EBs were first discovered by 
Ellerman in 1917 (jEUermanl 119171 ). The spectral profile of EBs consists of an absorption 
core at the Ha line center and asymmet r ic wide emissions in both red and blue wings of H a 
JSevernvill968l : lEngvold & Maltbvill968l : iKoval fc Severnv!ll97ol : lBruzeklll972l : lKitailll983h . 



The mean lifetime of EBs is es timated to be 10-20 minutes (IKurokawa et al.l 11982 
Nindos fc Zirinlll998l : iQiu et al.ll2000r). The typical s i ze of EBs is of the order of 1" and their 



shapes are mostly ellip tical (IKurokawa et al.l Il982l ; iNindos fc Zirinl Il998l : IGeorgoulis et al. 



2002: 



19721: 



Pariat et al 



Hu. Song, fc Li 



20071). The energy of EBs is estimated in the range 1 25 -10 28 ergs (iBruzek 



19951 : iHenoux et al1ll998l : IGeorgoulis eTaDbooi IFang et alibood ) 



It is re ported that EBs are associated wit h chromospheric upflows with a velocity about 
6-8 km s" 1 (IKurokawa et al.lll982l : lKitailll983l ). On the other hand, Georgoulis et al. (2002) 
found that almost all EBs were also accompanied by p hotospheric downflows with a yelocit y 



of 0.1-0.4 km s 1 . EBs are also related to Ha surges (jRustlll968l ; iMatsumoto et al. 1120081 ) 



which are thought to be evidence of magnetic reconnection in the low chromosphere. Some 
of chromospheric anemone jets observed with Ca II H filters are suggested to be associated 
to EBs (Shibata et al. 2007). 

It is important to study the relationship between EBs and magnetic fields, because EBs 
generally occur in areas of flux emergence and strong magnetic fields. In previous studies, it 
is reported that several EBs appear at the boundar i es of well-defined magnetic features, i.e., 
near m agnetic neutral lines (|Rov fc Leparskaslll973l : iDara et al.ll997l). near bipolar e merging; 
fluxes ( jZachariadis et al.lll987l ). or moving magnetic features (INindos &: Zirinl 119981 ) . In an 
emerging flux region (EFR), EB production is enhanced at the vicinity of sunspots, at areas 
near magnetic neutral lines and at the boundaries of the neighboring supergranular cells 
(IQiu et al.l l2000l : IGeorgoulis et al.ll2002l ). Georgoulis et al. (2002) suggested that EBs occur 
preferentially on separatrix or quasi-separatrix layers in the low chromosphere. 

Based on observational facts, several triggering models of EBs have been p roposed. Most 
generally accepted one is magnetic reconnect i on in the low chr o mosphere 



19981 : IGeorgoulis et al.l 12002 : IFang et al.l 120031 : IPariat et al.l 12004 : IFang et al 
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20061 ). Chen 



et al. (2001) confirmed that the magnetic reconnection in the lower solar atmosphere can 
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explain the main characteristics of EBs using a two-dimensional MHD simulation. On the 
other hand, Qiu et al. (2000) found that the majority of the Ha-UV(ultra violet continuum at 
1600A) well-correlated EBs were located at the boundaries of unipolar area, and argued that 
EBs located at unipolar magnetic areas could be triggered by another mechanism different 
from magnetic reconnection at low chromosphere. From topological analysis, Pariat et al. 
(2004) found that EBs and bald patches were linked by a hierarchy of elongated flux tubes 
showing aperiodic spatial undulations, whose wavelengths were typically above the threshold 
of the Parker instability. These findings led to a resistive emergence model of magnetic flux 
tubes: EBs are the signature of resistive emergence of undulatory flux tubes provoked by 
Parker instability. Isobe, Tripathi, & Archontis (2007) confirmed the validity of this model 
by a two-dimensional MHD simulation. 

Our motivation is to test magnetic reconnection model of EBs by analyzing time series 
data of magnetic field and velocity field in an EFR. In this paper, we use the high resolution 
data from spectro-polarimetric observation of NOAA 10917 taken with a vector magneto- 
graph at Domeless Solar Telescope (DST). Observational information and data reduction 
are presented in §2. The analysis of morphological characteristics from slit scan observation 
are described in §3. Temporal evolutions of magnetic field and velocity field around EBs 
observed in slit fixed observation are shown in §4. Finally in §5, we give a discussion and 
our conclusion about the triggering mechanism of EBs. 



2. OBSERVATION AND DATA REDUCTION 

We observed a (3 sunspot in an active region (AR) NOAA 10917 from October 21, 
2006 through October 22 for about 2.5 hours. Left panel in figure [T] shows a full sun image 
of October 21 taken with Solar Magnetic Activity Research Telescope (SMART) at Hida 
Observatory. With DST, we performed imaging observation of the region in Ha on band 
and off band, and Fe I spectroscopic observation with the vector magnetograph. The AR 
consisted of two main preceding and following spots, and several small spots (see right 
panel in figured]). In Ha line center images (figure QJc)), there were several dark filaments 
connecting the leading and the following one. The dark filaments showed red shift at their 
footpoints and blue shift at their tops, therefore they were rising flux tubes. Beneath the 
dark filaments, there appeared numerous EBs in Ha wing images (figure Q^a), (b)). The 
heliocentric coordinate of NOAA 10917 at 00:00 UT on October 22 was 5°S and 33°W. The 
AR was rather stable and in its developing phase during our observation period. No flare 
activity during our observation was reported by Space Environment Center. 
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2.1. Ha filtergram 

10 wavelengths set of Ha images, line center, ±0.3A, ±0.5A, ±0.8A, ± 1.2A, -5A, were 
obtained repeatedly with DST from October 21 22:50 UT through October 22 01:27 UT. 
In our analysis, we used Ha ±0.8A images for identification of EBs. The pixel size of Ha 
filtergram was 0".256, and the field of view (FOV) was 261" x 261". The time cadence of each 
wavelength frame was ~40 seconds. All the data were corrected by dark current subtraction 
and flat fielding. 

2.2. Vector magnetograph 

Along with Ha filtergrams, Fe I spectra were observed using a vector magnetograph 
(VMG). The spectral FOV includes both absorption lines of Fe I 6301.5A and Fe I 6302.5A. 
The spectral resolution is 0.0083A pixel" 1 at 6302.5A. The slit length was 140".8 and the 
width was 0".275. VMG consists of a rotating waveplate, three Wollaston prisms, and a 
high-dispersion spectrograph (Kiyohara et al. 2004). During one rotation of the waveplate, 
a pair of spectral images in orthogonally polarized states were taken at every 22.5° interval. 
It took ~30 seconds for one rotation of the waveplate. 16 pairs data were demodulated into 
the Stokes vector I, Q, U, V with the method described in Kiyohara et al. (2004). The 
calculated Stokes I, Q, U, V were inverted to parameters such as magnetic field intensity, 
magnetic inclination, magnetic azimuth, and line-of-sight velocity by ASP (Advanced Stokes 
Polarimeter) inversion code (see Skumanich & Lites (1987)). The ASP inversion technique 
adopts a hypothesis of a Milne-Eddington atmosphere, so that parameters are assumed to 
be constant along line-of-sight direction. Treatment of inversion was applied only in the 
area whose polarization degree was larger than 0.5%. Noise levels of the observed quantities, 
which were estimated from the temporal fluctuations at a relatively stable point near the 
center of the AR, were ±50 Gauss for magnetic field intensity, and ±3° for inclination and 
azimuth angle. 

We performed two kinds of observation mode using VMG. One was "slit scan" and 
the other was "slit fixed." In the slit scan mode, we scanned the AR by 2" intervals for 
46 positions, so that we could get one map with 92" x 140". 8 FOV. In the slit fixed mode, 
we tracked the AR manually aligning the slit fixed along the axis of the AR, where EB 
production was active. The exposure time and the time cadence were 1500 ms and ~30 
seconds for both observation modes. First, slit fixed observation was done from October 21 
23:03 UT to October 22 00:53 UT. Next, slit scan observation was done from October 22 
01:00 UT to 01:25 UT. 
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All spectral data were corrected by dark current subtraction and flat fielding. Slit fixed 
data were co-aligned with reference to the position of the main spots on continuum images. 
We performed box car smoothing of 5 frames (~2.5 minutes) on the slit fixed data to reduce 
high frequency noise. For time series data of Doppler velocity, we applied lowpass filter 
(<3 mHz) in order to eliminate the effect of photospheric five minutes oscillation. To get 
the physical insights into the magnetic configuration more easily, we corrected for geometric 
projection effect to the vector magnetic field. To be more precise, the magnetic inclination 
was converted to the tilt angle from local normal to the solar surface, and the azimuth was 
converted to the tilt angle measured from the slit direction on local horizontal plane. 

2.3. Identification of EBs 

As was stated in §1, EBs are characterized by their dark core and bright wing enhance- 
ments in Ha. To identify EBs, we calculated the contrast I c (x,t) = [I(x,t) — I (t)]/I (t) on 
Ha ±0.8A images. Here I(x, t) is the intensity of a pixel at position x taken at time t, and I 
is the mean intensity of the entire frame at time t. We adopted an empirical threshold value 
of I c =1.15 for EBs, that is, EBs are the regions where intensity is more than 15% brighter 
than the average. This value is consistent with earlier works, 5%— 30% in Georgoulis et al. 
(2002) and 14% in Pariat et al. (2007). We made sure that they did not show any brightness 
enhancements at the positions of EBs in Ha center images. 

In order to know the EFR magnetic conditions around EBs, we analyzed only such EBs 
that were adjacent to the VMG slit. As a result, 12 EBs were detected during slit scan 
observation, and 9 EBs during slit fixed observation. Figure [2] shows temporal series of Ha 
wing and line center images and the position of 9 EBs during slit fixed observation period. 
In figure [21 the appearances of the dark filaments and their velocity distributions change 
rapidly in time. This may be due to new emerging flux tubes (EFTs) which pressed the 
dark filament from below the photosphere. Among 9 EBs during slit fixed observation, 7 
EBs were observed from their birth to their death, another EB already existed when the 
observation started, and 1 EB did not yet decay when the observation stopped. 

3. RESULTS OF SLIT SCAN OBSERVATION 

Figure [3H6] show results obtained from slit scan observation. The slit scanning direction 
was in transversal, from left to right in each figure. The positions of 12 identified EBs are 
indicated by yellow or white diamonds. In this section, we mainly describe the morphological 
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characteristics of the AR and preferential locations of EBs. 



3.1. Magnetic field 

Figure [3] shows the distribution of magnetic field intensity. Small spots had large field 
intensity ~2000 Gauss (see figure [T] for reference). EBs were distributed at regions with 
field intensity 500-1000 Gauss. Field inclination angle shown in figure H] means the tilt angle 
from local normal to the solar surface, and the background gray scale in figure [5] shows 
vertical components of magnetic field, i.e., (field intensity) x cos (field inclination). After the 
correction for the projection effect, we solved 180° ambiguity of azimuthal angle with a 
hypothesis that the magnetic lines were similar to those of a dipole between the two main 
spots, because NOAA 10917 was a simple dipole in overview. Green arrows in figure [5] 
indicate horizontal vectors of magnetic lines whose lengths are proportional to log-scaled 
horizontal components of magnetic field at each point. As is seen in the background image 
of figure that EB production areas showed large spatial gradient of vertical field intensity. 
This suggests the possibility of rapid topological change of magnetic field in the area of EBs, 
probably due to new EFTs. 

We found that most EBs were distributed at regions whose inclination angles ~90°, 
that is, near magnetic neutral line. We found 8 out of 12 EBs near magnetic neutral lines 
(square and cross signs in figure 0H6]) , while 4 EBs in unipolar region. Here, we introduce 
two classificati ons of neutral lines , "top" and "dip." They are defined by the sign of the field 



line curvature (IPariat et al.ll2004t ). "Top" is a region of the photosphere (z = 0) where the 



field line curvature is negative, 

B z = and B • VB Z < (1) 

while "dip" is a region where the field line curvature is positive, 

B z = and B • VB Z > (2) 

"Top" is sometimes referred to Q loop, and "dip" is U loop. Then the 8 EBs found on neutral 
lines are classified into 5 "top" EBs and 3 "dip" EBs. In figure HHHl the cross symbols mean 
"top" and the squares mean "dip." 



3.2. Doppler velocity 



Figure [6] shows Doppler velocity field, or line-of-sight velocity field. Red regions have 
receding velocity from us and blue regions have approaching velocity to us. The effect of 
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solar rotation velocity is already subtracted. There was no significant velocity field with 
\v\ > 2 km s _1 around this AR. We assumed that the observed Doppler velocity consisted of 
large scale flow field and localized flow field. In large scale view, blue trend in the east side 
and red trend in the west side were considered to be horizontal flux segregation motion of the 
EFR or supergranular horizontal diverging flow (Bernasconi et al. 2002; Kozu et al. 2006; 
Niagara 2006). Because the target's position was in the west side, the horizontal segregation 
motion is observed as blue in the east side and red in the west side. Using SOHO/MDI full 
disk magnetogram images, we estimated the horizontal velocity around this AR based on the 
local correlation tracking method. The averaged horizontal velocity over our observational 
period was 0.16 km s _1 , diverging from the center to east and west direction. Then the 
effective component of horizontal velocity against the line-of-sight velocity was less than 0.1 
km s _1 . 

Let us look at the localized relation between EB positions and Doppler map in figure 
[6j The EBs were mostly located on the boundary of red/blue Doppler velocity. This result 
also proved the existence of EFTs. One example is shown in the areas enclosed with a 
green rectangle in figure H] and El There existed an Q shaped magnetic field line. Its top 
showed blue shift, or upward motion of the magnetic line and its footpoints showed red 
shift because of the plasma sliding down along the loops. Thus the observed local Doppler 
variation around the EBs can be interpreted as was sketched in figure [7J So the estimated 
amplitude of horizontal velocity (< 0.1 km s _1 from SOHO/MDI magnetogram images) and 
the results of earlier works (Lites et al. 1998; Kozu et al. 2006) allow us to suppose the 
localized line-of-sight velocities were predominantly contributed by the vertical motion. 



4. RESULTS OF SLIT FIXED OBSERVATION 

Next, we proceed to the results obtained from slit fixed observation. Figure [8] explains 
the coordinate system used in figure l9]fl~2l The left and right Ha images are respectively Ha 
— 0.8A images taken at the beginning of the observation and at the end. For example, in the 
left Ha panel, we can see a bright patch near the center of the slit, which is marked as EB1 in 
the middle panel. The middle panels show the temporal variation of physical parameters on 
the slit. The dashed and dash-dotted rectangles are indicators of EB positions. The dashed 
rectangles show EBs found by Ha — 0.8A images, and the dash-dotted ones show EBs found 
by Ha +0.8A images. As is demonstrated at upper right in each figure, the horizontal length 
of the rectangle indicates its lifetime, and the vertical width covers a zone of the EB position 
±3". 4. The mean lifetime of the 9 detected EBs was ~20 minutes, which is comparable to 
those of previous works. As is in the slit scan observation, 180° ambiguity of azimuthal angle 
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is resolved with a hypothesis that the magnetic lines were similar to those of a simple dipole. 

The slit fixed figures I9TTI21 appear to illustrate another important EB characteristic, i.e., 
that EBs appear to occur and recur at preferential locations. Because the slit's position was 
fixed along the axis of the AR, we can easily classify into three groups of EBs that faded and 
reappeared at nearly the same locations. These groups are (EB3, EB5, EB9), (EB6, EB7), 
and (EB1, EB4, EB8). This is a nice result which shows occurrence and recurrence of EBs 
with temporarily varying underlying magnetic conditions. 

4.1. Magnetic field intensity 

Figure [9] shows temporal evolution of magnetic field intensity on the slit. As we men- 
tioned in §3.1, EBs were distributed at regions with 500-1000 Gauss magnetic field intensity. 
Magnetic field intensity near EBs was evidently increasing, or at least changing in time. If 
a new EFT makes an appearance from beneath the photosphere, it presses the pre-existing 
magnetic field and magnetic intensity gradually increases. Therefore strong magnetic inten- 
sity area can be interpreted as an EFT production site. Except for EB6, EB7, and EB9, the 
other 6 EBs showed signatures of increasing field intensity. 

4.2. Magnetic inclination 

In figure [TDJ temporal evolution of magnetic field inclination with respect to local nor- 
mal (+z) is shown. We made a more detailed picture of temporal variation of inclination 
angle around each EB in figure [13j In figure Q3J arrows indicate the magnetic vectors 
(y,z) = (cos(inclination), sin(inclination)) with unit length. As is the same in §3.1, almost all 
the EBs were distributed near magnetic neutral lines. 

According to the definition explained in §3.1, EB1, EB2, EB4, EB6, EB7 and EB8 
belonged to "dip." No "top" EB was found. And EB3, EB5, and EB9 belonged to neither 
of them but on unipolar region. Although EB3 and EB5 did not lie in magnetic neutral 
line but stay in unipolar region, their distributions of inclination angles showed apparent dip 
configuration, which is clearly seen in figure [T3J We shall call this kind of configuration as 
"local dip" from now on. 

At the locations of EB1 and EB3, we found temporal formation of magnetic dip struc- 
ture. At first, the distribution of magnetic field was uniform along the slit. As the time 
passed, inclination of field lines evolved such that the magnetic topology took the dip struc- 
ture in the photosphere. It is interesting that Ha wing brightenings already started while the 
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magnetic field lines were still in uniform state in the photosphere, which may indicate the 
magnetic dips were already formed in the upper layers than the Fe I 6302. 5A line formation 
heights and started to produce an EB by magnetic reconnection mechanism. 

In the central part of the active region, there was a clear undulatory pattern, which 
is seen as periodically blue and red regions along the slit in figure [10J All the EBs were 
distributed in this undulatory field line area. This is supporting evidence for the resistive 
emergence model of Pariat et al. (2004). 

4.3. Magnetic azimuth 

Figure [TT1 shows temporal evolution of magnetic field azimuth on local horizontal plane. 
Azimuthal angle is the tilt angle measured from the slit direction (+y). Magnetic lines 
in orange colored region are inclined to right (+x) direction and magnetic lines in green 
colored region are inclined to left (— x) direction. Azimuthal angles between the two main 
spots mostly took value ranging from —20° to +20°. Like inclination angle, azimuthal angle 
showed a clear undulatory pattern. Two plots of inclination and azimuth angle along the 
slit are shown in figure [TH The four cross signs are indicators of EB generation sites. It is 
clearly seen that the magnetic field around EBs were undulatory, and that EBs were located 
in the region where the inclination of magnetic lines was ~90° and their azimuth angles 
changed their signs. The characteristic wavelength of the undulation pattern exceeded 5", 
that is, longer than 4000km. This result tells us an important fact, because 4000 km is 
the characteristic wavelength of Parker instability at the photospheric layer. As long as we 
know, this is the first to report an undulatory pattern in azimuthal angles around EBs. 

4.4. Emerging flux tubes 

In the observed AR, there were signatures of many EFTs in the photosphere. The pref- 
erential position of EBs showed an undulatory pattern in their magnetic field components, 
probably due to new EFTs triggered by Parker instability. We identified three small EFTs 
during the slit fixed observation. Figure [TBI shows the positions of three EFTs and EBs. The 
identification of an EFT was done by examining the distribution of magnetic field intensity, 
inclination, and Doppler velocity. Figure [16] shows three kinds of plots, magnetic field inten- 
sity, inclination and Doppler velocity along the black solid line in figure l9lfT2l The field line 
showed Q shape in its inclination, the field intensity was increasing, and Doppler velocity 
changed its sign across the neutral line. Blue shifted area was found near the top of the 
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Q loop field lines. Blue shifted area can be interpreted as vertically upward motion as was 
suggested in §3.2. At both sides of the footpoints, there existed red shifted, i.e., downward 
moving areas. We think these downflows were due to the motion of plasma dragged from 
the top of the Q loop along field lines. Then we conclude this area as an EFT, and call this 
EFT1. EFT1 was located near EB1 and started emergence from about 23:10 UT and had 
not yet decayed at the end of our observation. Second one (EFT2) was located near EB2. 
It started emergence before the observation started and decayed at about 00:00 UT. EFT2 
also showed magnetic intensity enhancement and Doppler blue shift ~0.1 km s -1 at the top, 
and Doppler red shift ~1.0 km s _1 at its footpoints. Third one (EFT3) made an appear- 
ance near EB3. It started emergence at about 23:20 and had not yet decayed at the end of 
the observation. EFT3 induced magnetic field intensity enhancement at its early phase and 
Doppler blue shifted velocity ~1.0 km s _1 . However unlike EFT1 and EFT2, EFT3 did not 
have red shifted areas at its footpoints. 



5. DISCUSSION 

5.1. Schematic models 

With slit scan observation, we found that 8 out of 12 identified EB production areas 
had large spatial gradient of vertical field components. This indicates the existence of new 
EFTs. Doppler velocity map in figure [6] shows that almost all the EBs had downward motion 
adjacent to upflow, which is also a signature of new EFTs. Thus we confirmed that EBs 
have a strong relation to the emergences of EFTs with slit fixed observation. As for the 
residual 4 EBs found on unipolar region, we also try to interpret by the emergence of EFTs 
and magnetic reconnection, as will be discussed below. 

During the slit fixed observation period, we identified three EFTs (figure [To^ . Each 
EFT invoked a gradual enhancement of magnetic field intensity. In regions where magnetic 
fluxes emerge continuously, new EFTs compress the pre-existing fluxes and the flux density 
increases. If oppositely directed or at least sheared magnetic lines are compressed and 
approaches each other, a magnetic reconnection occurs. It is widely accepted that EBs are a 
signature of magnetic reconnection at the low chromosphere. Therefore, it is plausible that 
magnetic field intensity increases at the location of an EFT and an EB. 

In figure [T7| we show schematically how our data are interpreted by magnetic reconnec- 
tion. Model 1,2, and 3 show EB mechanisms with different magnetic configurations: Model 
1 for footpoint of an EFT, model 2 for top of an EFT, and model 3 for unipolar region. 
According to the definition of magnetic lines in §3.1, a "dip" EB is explained by model 1, a 
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"top" EB is explained by model 2, and a unipolar by model 2 or model 3. Due to the coarse 
scan step and lack of temporal information, we can not clarify that 4 unipolar EBs found on 
the slit scan observation were located at the top of an EFT or not, and so 4 unipolar EBs 
can be explained either of model 2 or model 3. In summary, the identified EBs are classified 
as follows; (EB1, EB2, EB4, EB6, EB7, EB8) and 3 "dip" EBs of slit scan belonged to model 
I. (EB3, EB5, EB9) and 5 "top" EBs of slit scan belonged to model 2. Finally, 4 unipolar 
EBs of slit scan belonged to model 2 or model 3. 

In model 1, new EFTs appear from below the pre-existed uniform magnetic field. There 
will be produced an anti-parallel layer at the low chromosphere, just above the photospheric 
"dip" region, and reconnection occurs at X point. When reconnection occurs, there should 
exist converging flow across anti-parallel field lines. This speculative flow is shown in figure 
[17] with broken line arrows. We could not detect this speculative horizontal flow in our obser- 
vation, probably due to their small amplitude of velocity, compactness of spatial extension 
of the velocity field, and etc. 

Model 2 explains the triggering mechanism of an EB on top of an EFT. Model 2a is the 
case for local dip configuration, which accounts for EB3 and EB5. Like in model 1, horizontal 
converging flow may exist across the local dip position (broken line arrows). Model 2b is 
the case for no local dip configuration. If there is a shear between the neighboring magnetic 
lines, magnetic reconnection can occur at that point. In fact, Linton & Antiochos (2005) 
reported that reconnection can occur even at small contact angles between the interacting 
magnetic field lines. In model 2b, only upward motion will be observed. This may account 
for EB9. The 6 "top" EBs and 4 unipolar EBs on slit scan observation can be understood 
by both model 2a and model 2b. 

In model 3, there is no dips of magnetic field lines. If there is a shear between the 
neighboring magnetic lines, magnetic reconnection occurs at that point. This may account 
for unipolar EBs on slit scan observation. We are not sure how velocity field is distributed 
in model 3. 

As for EB1 and EB3, inclination angles were at first uniform and gradually made dip 
configurations. These two EBs appeared at the early phase of EFTs. The dip configuration 
may be observed only after the EFT has developed enough, because of the spatial resolution 
or the formation height. So even if we have a snapshot of magnetic inclination angle around 
EBs that do not show dip configuration, we can not exclude the possibility that this EB is 
not related to dip areas. 

Pariat et al. (2007) reported that EBs are located in regions named as bald patches, 
where the filed lines have geometrical form of U shape. This result is somehow consistent 
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with our model 1 and model 2a. We extended their work by observing temporal evolution 
of magnetic configuration of an EFT. In our data, more than half of the observed EBs can 
be accounted for by model 1 or model 2a. However, there were a non-negligible number of 
the observed EBs that were located at the top of the Q loop or unipolar region. Therefore 
model 2b and model 3, that is, EB triggering without dip configuration are also the effective 
mechanisms for EBs. 



5.2. Evidence of Parker instability 

We found an undulatory pattern whose characteristic wavelength was longer than 4000 
km on both inclination and azimuth data (figure HI]) . In Pariat et al. (2004), the mechanism 
of an EB is ma gnetic reconnection of undulatory flux tubes emerging by the Parker instability 



( jParker I Il966l ). If the undulatory magnetic fluxes are emerged by the Parker instability, 
magnetic field lines become undulatory in its vertical components. Their wavelength A 
should satisfy 

A > 4ttH ~ 2000 km (3) 

Here H is the pressure scale height, and take value of 150-200 km at the photosphere. This is 
called a resistive emergence model. The undulatory pattern in magnetic field inclination is 
direct evidence for this model. In addition, the undulatory pattern in magnetic field azimuth 
also prove Parker instability by the following reason. As the magnetic fluxes that emerge by 
the Parker instability go up, the atmospheric pressure decreases as height and the emerging 
flux expands horizontally (Shibata et al. 1989). But the footpoints anchored in photosphere 
are still compressed by surroundings. When observed this EFT from the top, the horizontal 
components of magnetic field line shows wavy undulations with a characteristic wavelength 
of the Parker instability. This is the first to report an undulation of horizontal components 
of magnetic field, which is supporting evidence for Parker instability. 
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Fig. I. — Left: Full sun image in Ha line center on October 21, 2006 taken with SMART. 
NOAA 10917 is pointed with an arrow. The square box indicates the field of view of right 
panels. Right: Ha images of NOAA 10917 taken with DST. Four images are taken nearly 
simultaneously at October 22 00:59 UT. (a)Hct -0.8A (b)Hct +0.8A (c)Hct line center (d)Ha 
— 5. OA. Three emerging flux tubes (EFTs) which show downward motions at their footpoints 
and upward ones at their tops are encircled with broken lines. 
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Fig. 2. — Ha wing and center images taken during the slit fixed observation. The upper 6 
panels are Ha +0.8A images, in which downward moving regions appear dark. The middle 6 
panels are Ha center images. Lower 6 panels are Ha -0.8A images, in which upward moving 
regions appear dark. The position of each identified EBs are pointed with arrows. 




Fig. 3. — Spatial distribution of magnetic field intensity (slit scan). White diamonds show 
the positions of observed EBs. 
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Fig. 4. — Field inclination with respect to local normal (slit scan). Yellow diamonds show 
the positions of observed EBs. Cross signs indicate EBs at the "top", and square signs 
indicate EBs at the "dip". The magnetic configuration in green rectangle area is given in 
figure [7J 
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Fig. 5. — Vertical component of magnetic field intensity (background gray scale) and horizon- 
tal magnetic field vectors (slit scan). The length of each arrow is proportional to log-scaled 
horizontal components of magnetic field. Yellow diamonds show the positions of observed 
EBs. Cross signs indicate EBs at the "top", and square signs indicate EBs at the "dip". 
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Fig. 6. — Line-of-sight velocity (slit scan). Negative velocity means vertically upward or hor- 
izontally eastward motion, and positive velocity means vertically downward or horizontally 
westward motion. Yellow diamonds show the positions of observed EBs. Cross signs indicate 
EBs at the "top," and square signs indicate EBs at the "dip". The magnetic configuration 
in green rectangle area is given in figure [71 
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Fig. 7. — Plausible configuration of magnetic field lines in the area enclosed with a green 
rectangle in figure H] and figure [6j Blue and red arrows indicate the direction of plasma mass 
flow. 




Fig. 8. — Coordinate system used in slit fixed observation results. 
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Fig. 9. — Magnetic field intensity variation with time (slit fixed). Times are shown in 
UT. The rectangles are indicators of EB positions. The horizontal length of the rectangle 
indicates its lifetime, and the vertical width covers a zone of the EB position ±3". 4. Two 
Ha -0.8 A images are given for references at the beginning and the end of the observation. 
Plot along the black solid line is shown in figure fToTa). 
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Fig. 10. — Magnetic inclination variation with time (slit fixed). The others are the same as 
in Fig. [9j Plot along the black solid line is shown in figure [ToT b). 



- 25 




Fig. 12. — Doppler velocity variation with time (slit fixed). Negative means blue shift, and 
positive means red shift. The others are the same as in Fig. [91 Plot along the black solid 
line is shown in figure [ToT c). 
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Fig. 13. — Temporal evolution of magnetic field inclination around each EB. Time for each 
frame is shown at upper left in UT. Circle symbol at the center of each figure indicates the 
position of the EB. 
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Fig. 14. — (a)Magnetic inclination and (b)azimuth along slit at 23:33:32 (slit fixed). Four 
cross signs are the indicators of the position of EB occurrence. Position i (EB3, EB5, EB9). 
Position ii (EB6, EB7). Position iii (EB1, EB4, EB8). Position iv (EB2). 




Fig. 15. — Position of three EFTs and EBs. The ground figure is the same as the left Ho; 
image shown in figure IMT21 
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Fig. 16. — Three kinds of plots of EFT1 located along a black solid line in fi gure l9l [TUl and [T2l 
(a)magnetic field intensity (b)magnetic inclination (c)Doppler velocity (Negative is upflow 
and positive is downflow.) 
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Fig. 17. — Three different mechanisms of EB triggering, (model 1) EB mechanism at the 
footpoints of EFTs. (model 2) EB mechanism at the top of EFTs. (model 3) EB mechanism 
at unipolar region. The layer at the footpoint of vectors indicates the formation height of Fe 
I 6302. 5A. X points are reconnection points. Lower rectangles with minus symbol mean that 
inclination angle is larger than 90° and rectangles with plus symbol mean that inclination 
angle is less than 90°. Thick arrows show Doppler velocities; blue arrows for upflow and red 
arrows for downflow. Broken line arrows written in model 1 and model 2a means horizontal 
converging flow that may exist. 



